Crossover from the pair contact process with diffusion to directed percolation 
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Crossover behaviors from the pair contact process with diffusion (PCPD) and the driven PCPD 
(DPCPD) to the directed percolation (DP) are studied in one dimension by introducing a single 
particle annihilation and/or branching dynamics. The crossover exponents are estimated numeri- 
cally as 1/0 ~ 0.58 ±0.03 for the PCPD and 1/0 ~ 0.49 ±0.02 for the DPCPD. Nontriviality of the 
PCPD crossover exponent strongly supports the non-DP nature of the PCPD critical scaling, which 
is further evidenced by the anomalous critical amplitude scaling near the PCPD point. In addition, 
we find that the DPCPD crossover is consistent with the mean field prediction of the tricritical DP 
class as expected. 

PACS numbers: 64.60.Ht,05.70.Ln,89.75.Da 



The absorbing phase transition (APT) has emerged 
during the last few decades as a prototype of nonequi- 
librium critical phenomena. The APT is a transition 
from an active phase into an inactive (absorbing) phase 
which is composed of absorbing (trapped) configurational 
states where the system cannot leave by the prescribed 
dynamic rules. As in the equilibrium critical phenom- 
ena, a few scaling exponents characterize and classify 
critical behaviors into universality classes 0, Q. Ro- 
bustness of the directed percolation (DP) universality 
class with respect to the microscopic details led to the 
"DP conjecture" 0, that a model should belong to 
the DP class if it has a unique absorbing state without 
additional symmetry, conservation laws, quenched dis- 
order, and long-range interactions. Although significant 
progress has been achieved to date, the full understand- 
ing on the main features of the APT universality classes 
is still far from complete. 

To achieve such a goal, it is crucial to identify the uni- 
versality class of the pair contact process with diffusion 
(PCPD) HSI1, 

which has been the most controver- 
sial topic these days, to our knowledge. The PCPD is a 
variant of the pair contact process Q by allowing single- 
particle diffusion in addition to binary fission and anni- 
hilation dynamics (2 A — > 3 A and 2 A — > 0). Even with 
combined efforts of extensive and highly equipped numer- 
ical and analytical studies, its one-dimensional version 
has defied a consensus as yet @. Rather, various possi- 
ble scenarios have been suggested with the DP class in 
the center of this controversy. In one scenario, the PCPD 
should belong to a universality class distinct from the DP 
with a unique set of critical exponents 0, [ToL ITl| or con- 
tinuously varying exponents due to the marginal pertur- 
bation 12] . The other states that the PCPD will even- 
tually flow to the DP fixed point after a huge crossover 
time 0, 0, 0| . More scenarios can be found in Ref . @ . 
In high dimensions, it is unquestionably clear that the 
PCPD and the DP exhibit different critical scaling. 

The strong corrections to scaling are the main obsta- 
cle for the numerical study to determine the universality 
class of the PCPD. To make matters worse, the field the- 
ory with a single component order parameter is shown to 



be not viable |l6j | . It implies that the proper field theory, 
if it exists, needs (at least) two independent order param- 
eters, which is also independently proposed by us from 
the numerical study of the driven PCPD (DPCPD) [To| . 

In this Rapid Communication, we propose an approach 
to settle down or at least moderate the PCPD contro- 
versy by scrutinizing closely the crossover behavior from 
the PCPD to the DP class. We allow unary branching 
and/or annihilation dynamics (A — > 2A and A — ► 0) in 
the PCPD model, which should lead to the well-known 
DP critical dynamics as soon as it is introduced. The 
nontrivial nature of the DP critical line approaching the 
PCPD point as well as anomalous crossover scaling be- 
havior would serve as an indicative of different scaling at 
the PCPD point, distinct from the DP. 

Although the accurate estimation of the critical expo- 
nent values for the PCPD is severely plagued by correc- 
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FIG. 1: (Color online) Log-log plot of the particle density 
decay for the PCPD at p = 0.133 516. After t ~ 10 s , the 
density saturates, which signals that the system is in the ac- 
tive phase. Inset: Semilogarithmic plot of p(t)t s versus t at 
p = 0.133 522 with S = 0.19. After t ~ 10 7 , the curve veers 
down, which signals that the system is in the absorbing phase. 
This result is insensitive to the change of the S value by up 
to 15%. 



2 



tions to scaling, the critical point itself can be located 
with relatively high precision. In Fig. ^ the density 
decay near criticality is shown for the PCPD model of 
which detailed dynamics is given in Eq. Q with w = 0. 
The system size at a tuning parameter p = 0.133 516 
(0.133 522) is 2 18 (2 20 ) and the number of independent 
samples is 32 (104). All samples carry at least one neigh- 
boring pair until the end of simulations, t = 10 9 (10 8 ), 
which guarantees that finite size effects have not emerged 
as yet. From Fig. ^ the critical point is estimated as 
p c = 0.133 519(3) with the number in the parentheses 
as the uncertainty of the last digit. It should be em- 
phasized that our estimation of the critical point for the 
PCPD does not resort to any prescribed critical exponent 
value. 

In the generalized crossover model including unary 
branching and/or annihilation processes, it is relatively 
easier to locate the critical points accurately by utiliz- 
ing the known DP critical exponent values. With the 
annihilation A — > and/or the branching A — > 2 A, the 
density decay in the absorbing phase should be exponen- 
tial in time at least in one dimension, which nullifies the 
long-term memory effect of the isolated particles 0] an d 
the simple unplagued DP scaling should be anticipated. 
The accurate information on the critical line as a func- 
tion of the unary process rate w allows us to estimate 
the crossover exponent </>, which describes the crossover 
scaling, if it exists, from the PCPD to the DP. 

The detailed evolution rules for the generalized 
crossover model are summarized using stoichiometric no- 
tations as 



AA — > 
9AAj 




3 with rate A, (la) 
* AAA, with rate a/2, (lb) 

9 A with rate D(l-w)/2 

with rate w(\ — q)/2 , (lc) 
with rate wq/2 

with rate D(l - w)/2 
with rate w(l -q)/2 , (Id) 
with rate wq/2 



where A (0) stands for a hard core particle (vacancy) 
and < q < 1. The periodic boundary conditions are 
employed on a one-dimensional lattice with size L. The 
case of w = 0, D = 1, A = p, and a = 1 — p with 
a tuning parameter p corresponds to the PCPD model 
studied before °f which the critical point is accu- 
rately located through Fig. ^ I n the case of nonzero w, 
the critical point is located by observing the flatness of 
p(t)t Sop as a function of t with the DP critical exponent 
(5dp — 0.1595 |17|. In the following, we always set D = 1 
and A = 1 — a = p. 

The simulation algorithm to mimic the dynamics in 
Eq. (JTJ is as follows: First, choose a particle at ran- 
dom. Then, choose one of its nearest neighboring sites 
randomly as a target site. If the target site is vacant, the 
chosen particle is annihilated with probability w(l — q), 



TABLE I: Critical point values p c o(w) and p c i(w) for various 
w's for q — and q = 1, respectively. The numbers in the 
parentheses indicate the uncertainty of the last digits. 



w 


Pco(w) 


Pel (w) 





0.133 519(3) 


0.133 519(3) 


1(T 5 


0.133 172(3) 


0.134 085(5) 


5 x 1(T 5 


0.132 581(2) 




1(T 4 


0.132 080(3) 


0.135 81(1) 


2 x 1(T 4 


0.131 325(5) 


0.136 98(1) 


3 x 1(T 4 


0.130 715(5) 


0.137 92(1) 


4 x 1(T 4 


0.130 185(5) 


0.138 73(1) 


5 x 1(T 4 


0.129 708(3) 


0.139 46(1) 



branches a particle at the target site with probability wq, 
or hops to that site with probability 1 — w. The random 
selection of one of nearest neighbors amounts to the fac- 
tor | of the transition rates in Eqs. lfTc|) and (|ld|) . In 
the case where the target site is already occupied, both 
particles are annihilated with probability p or one ex- 
tra particle is generated at a randomly chosen nearest 
neighbor site of the pair with probability 1—p when that 
site is vacant. If the chosen site is already occupied, this 
branching attempt is rejected. After this update, the 
time increases by 1/N t , where N t is the total number of 
particles at time t just before the update. 

By observing how p(t)t Sop behaves in the asymptotic 
regime, the critical point values p c for nonzero w's are 
estimated. If the system is in the active (absorbing) 
phase, those curves should veer up (down). Only at crit- 
icality, a flat straight line can be observed. Initially, all 
sites are occupied. The simulations were performed up to 



t = 10 —10 and no finite size effects have been observed 




FIG. 2: (Color online) Semilogarithmic plots of p(t)t SoF 
as a function of t for the generalized crossover model at 
criticality. The upper (lower) three curves correspond to 
w = 5 x 10~ 4 , 10" 4 , and 10" 5 from top (left) to bottom 
(right) at q = (1). The middle two curves are redrawn us- 
ing the same data as in Fig. 0for the PCPD (w = 0) near 
the critical point. 
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FIG. 3: (Color online) Critical points p c {w) for various w's 
and their fitting functions. The lower (upper) curve corre- 
sponds to q — (1). The error bar is smaller than the symbol 
size. The curves represent the fitting results of the critical 
points with the function A c (w) = aw 1 ^ . The estimated 
crossover exponent is l/<j> = 0.58 ± 0.03. Inset: Log- log plot 
of A c (w)| vs w. The slope of the straight line is l/(j>. 
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FIG. 4: (Color online) Phase boundary of the DPCPD 
crossover model with llc'll . Just as in Fig. |3] the crossover 
exponent is estimated by the least-square fitting as 1/0 — 
0.49 ± 0.02 which is consistent with the mean field crossover 
exponent [Tsj| . Inset: Log-log plot of |A c (w)| vs w. 



the dynamics (fTcf) is replaced by 



with system size L — 2 18 — 2 20 . The numerical estimates 
for the critical point values p c o(w) at q = and p c i(w) at 
q = 1 arc tabulated in Table|U Figure|2]shows the critical 
density decay for various values of w at q = and 1. The 
flatness of the critical curves can be seen quite early even 
for very small w, while the PCPD lines (w ~ 0) do not 
show any flatness within our observation time. It may be 
understood either by a tremendously long dynamic cor- 
rection to the DP scaling 0] or by a non-DP scaling for 
the PCPD. In cither case, it is evident that there exists 
a diverging crossover time scale as w — > 0. 

The crossover scaling function near the PCPD critical 
point should take the form 

p(w, A; t) = t- s F{Aw~ 1 /' t ',tw^), (2) 

where = v\\ /0 with the crossover exponent 0, v\\ , and <5 
are the critical exponents of the PCPD, and A = p— p c {0) 
with p c (0) being the PCPD critical point value. The 
crossover time diverges as w — > such that r cross ~ . 

The crossover exponent can be calculated without 
the knowledge of the values of v\\ and 5 by studying how 
the DP critical line p c {w) approaches the PCPD criti- 
cal point. This line should be one of the renormaliza- 
tion group flow lines such as A c {w) ~ aw 1 ^ 0, where 
A c (w;) = p c {w)—p c {0) and a is a (nonuniversal) constant. 
In Fig. the crossover exponent is estimated using the 
least-square fitting, which turns out to have a nontrivial 
and universal value of 1/0 = 0.58 ± 0.03. This nontriv- 
ial nature of the critical line (0^1) signals strongly a 
possible non-DP scaling at the PCPD point. 

The similar analysis can be applied to the DPCPD 
with diffusion bias by modifying the dynamics in Eqs. 
Q in such a way that the dynamics (|ld|) is absent and 



A0 i ® A With rate ^ ~ W ^ {Id) 
[ 00 with rate w 

With w = 0, this model is the DPCPD studied in Rcf. 
[To). Since the diffusion bias should not change the DP 
scaling ^(j> this model with nonzero w is still expected 
to belong to the DP class. Figure 0] shows how the 
DPCPD crosses over to the DP as w increases and finds 
the crossover exponent 1/0 ~ 0.49 ± 0.02, which is con- 
sistent with the mean field prediction of the tricritical 
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FIG. 5: (Color online) Density-pair density ratio R{t) for 
various w's with q — 0.7 in the semilogarithmic scale. The 
critical points are located at p c = 0.134 120(5), 0.133 761(2), 
0.133 616(2), and 0.133 563(2) for w = 10~ 4 , 2 x 10~ 5 , 4 x 
10 -6 , and 10 -6 , respectively. The longest curve corresponds 
to the PCPD at p = 0.133 516. The saturating values are 
almost identical for all w's including the PCPD. Inset: Plots 
of R{t) for w's in Fig. [5] The upper (lower) three curves 
corresponds to q — 1 (0). For both cases, the ratio of two 
critical amplitudes approaches to the PCPD value as w be- 
comes smaller. 
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FIG. 6: (Color online) Log-log plots of p(t)t s ° F vs t for vari- 
ous w's as in Fig. I^Jwith q — 0.7 at criticality. Inset: Log-log 
plot of the critical amplitudes A(w) vs w. The slope of the 
straight line is 0.03. 

DP class [Tflj l and, in turn, confirms again that the upper 
critical dimension of the DPCPD is 1 poj ]. 

To get a more concrete evidence for the non-DP nature 
of the PCPD scaling, we go back to the scaling function 
in Eq. J2J) which can be rewritten at criticality as 

t Sop p(w,A c ;t) ^t Svp - s F(a,w^t) = w x g(w^t), (3) 

where \ = — fep) and Q(x) — x Sup ~ s F(a, x). Since 
Q{x) should be a finite constant as x —> oo to guarantee 
the DP scaling for nonzero w, the critical decay ampli- 
tude is proportional to w x for sufficiently small w where 
the crossover scaling is valid. 

If the PCPD belongs to the DP class (S = <5 D p), the 
amplitude converges to a nonzero value as w — » 0. Oth- 
erwise, it decreases algebraically with \ — 0.03, using the 
present best estimates for the PCPD exponent values as 
6 ~ 0.19 and v\\ ~ 1.85 HHQ1. In Fig 'El there 

are ap- 
parently two very different convergent behaviors for q = 



and 1. Naive estimates lead to x — 0.08 for q = and 
X — 0.00 for q = 1, both of which do not simultaneously 
fit into either the DP or the non-DP scenario. 

We find that this perplexing result originates from the 
narrowness of the crossover scaling region near the PCPD 
point. To get a hint for the valid scaling regime, we ob- 
serve the pair density p p as well as the particle density p 
in simulations. Both quantities should scale in the same 
way and their critical amplitudes also scale similarly as 
~ w x . Hence, the ratio R(t) = p{t)/p p {t) should be in- 
dependent of w for the valid crossover scaling regime in 
the asymptotic limit. The inset of Fig. [S] shows that 
both cases at q — and 1 have not entered the scal- 
ing regime as yet up to w = 10 -5 , which implies that 
the naive estimates for x should be significantly influ- 
enced by the corrections to scaling. Fortunately, we find 
a reasonably good scaling regime for not so small w at 
q = 0.7. Figurc[S]shows that the scaling regime is reached 
already for w = 1CP 4 . The inset of Fig. |Tj] shows how 
the critical amplitude of the density, say A(w), behaves 
in the crossover scaling regime, from which we estimate 
X = 0.03(1) that is consistent with the PCPD estimate. 
This result strongly supports again the non-DP nature 
of the PCPD scaling. The crossover exponent estimated 
for the case of q = 0.7 is 1/0 ~ 0.57, which is compatible 
with the previous estimation within error. 

To conclude, we studied the crossover behavior from 
the pair contact process with diffusion and the driven 
pair contact process with diffusion to the directed per- 
colation. We found that the crossover for the PCPD to 
the DP is nontrivial and the critical amplitude scaling 
is anomalous, which strongly supports the non-DP na- 
ture of the PCPD scaling. The crossover exponent for 
the DPCPD takes the mean field value of the tricritical 
DP (TDP) [lj|. This implies that the two-dimensional 
PCPD might have connections to the TDP that is under 
our current investigation. 
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